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摘要 : 高 空 核 爆炸 碎片 云 


了 


符 放 出 的 缓 发 Y 射线 在 高 空 非 均匀 大 气 的 输 运 过 程 中 电离 大 气 ， 使 得 电离 层 中 电 
子 数 密度 剧 增 ， 进 而 影响 途径 电离 层 的 无 线 电 通信 链 路 。 为 了 准确 描述 随时 空 演 化 的 组 发 Y 射线 源 电离 大 


气 过 程 ， 首 先 建立 高 空 核 爆 炸 碎片 运动 演化 的 流体 力学 模型 ， 然 后 根据 碎片 去 运动 参数 建立 组 发 Y 射线 分 


层 等 效 体 源 模型 ， 最 后 采用 蒙特 卡 罗 方 法 模拟 组 发 ?出 线 在 高 空 大 气 中 输 运 并 电离 大 气 的 过 程 。 针 对 碎片 云 
形状 随时 间 不 断 演 化 特性 ， 采 用 辐射 源 分 层 抽样 的 方法 得 到 缓 发 Y 射线 的 初始 位 置 。 针 对 大 气 密度 随 高 度 


指数 衰减 的 特性 ， 采 用 质量 厚度 抽样 方法 简化 计算 模型 。 结 果 表 明 : 碎片 云 形状 对 组 发 Y 射 线 电离 大 气 范 


围 和 强度 存在 明显 影响 。 兆 吨 级 高 空 核 爆 炸 缓 发 Y 射 线 电离 大 气 范 围 可 达 王 公里 以 上 。 当 爆 高 不 变 当 量 增 


加 时 ， 碎 片 云 高 度 和 水 平 半 径 增 大 ， 绥 发 Y 射线 电离 大 气 范围 和 强度 有 所 增 大 。 当 当量 不 变 爆 高 增加 时 ， 
碎片 云 高 度 和 水 平 半 径 增 大 ， 绥 发 Y 射 线 电离 大 气 范围 增 大 ， 但 强度 有 所 降低 。 
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Abstract: The delayed y-rays released from the debris of a high-altitude nuclear explosion ionize 
the atmosphere during transport in the high-altitude non-uniform atmosphere, causing a dramatic 
increase of the electron number density in the lonosphere, which affects the radio communication 
links through the ionosphere. In order to accurately describe the delayed y-ray sources which 
evolves in time and space, we firstly establish a hydrodynamic model of the evolution of debris 
motion of a high-altitude nuclear explosion, then establish a hierarchical equivalent model of 
delayed y-ray sources according to the debris parameters, and finally apply the Monte Carlo 


method to simulate the ionization effect of the delayed y-rays in the high-altitude non-uniform 


atmosphere. Because the debris shape continuously evolves with time, a _ method of stratified 
sampling of the radiation source 1s built to obtained the mitial positions of the delayed y-rays. Due 
to the atmosphere density decaying exponentially with height, the mass thickness sampling 
method 1s used to simplify the calculation model. The results show that the lonlzation intensity 
and range of delayed y-rays are significant influenced by the shape of the debris. The ionization 
range of delayed y-rays in a megaton level high-altitude nuclear explosion can reach more than a 
thousand kilometers. When the equivalent increases with a constant explosion height, the debris 
height and horizontal radius increases, the ionization range and intensity of the atmosphere caused 
by the delayed y-rays increase. When the burst height increases with a constant equivalent, the 
debris height and horizontal radius increases, the ionization range caused by the delayed y-rays 
increases, while the ionization intensity decreases. 
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高 空 核 爆 炸 中 产生 的 各 种 辐射 源 可 以 电离 高 层 大 气 ， 大 范围 长 时 间 增 加 电离 层 D、 王 层 
电子 数 密度 ， 破 坏 正常 电离 层 结构 53， 从 而 对 穿越 电离 层 的 无 线 电 波 造成 吸收 和 折射 等 多 
种 效应 ， 对 无 线 电 通 信 、 和 雷达 等 系统 带 来 影响 上 0]。 核 爆炸 所 释放 出 来 的 电离 辐射 按照 时 间 
的 先后 可 分 为 两 类 : 一 类 是 瞬 发 辐射 (包括 义 射线 、 中 子 、y 射 线 等 ) ， 是 在 核 爆炸 发 生 的 
瞬间 释放 出 来 的 ， 占 据 了 核 爆 炸 的 大 部 分 能 量 ， 因 此 可 以 在 大 气 层 中 形成 很 高 的 电离 。 另 一 
类 是 绥 友 辐射 (主要 是 y 冉 线 、B 粒 子 等 ) ， 是 由 核 爆 炸 裂 变 产 物 俯 片 云 氏 9 在 很 长 一 段 时 间 
内 逐渐 释放 出 来 的 。 相 比 于 瞬 发 辐 冉 源 , 组 发 辐射 源 是 随时 间 持 续 变 化 的 具有 明显 空间 分 布 
的 体 源 ， 且 持续 时 间 长 、 履 盖 范 围 大 。 

碎片 云 时 空 分 布 是 确定 绥 发 辐射 源 的 关键 基础 。 当 爆炸 高 度 和 当量 不 同时 ,碎片 云 形状 
可 以 由 球形 演化 为 扁 椭 球形 或 者 倒 梨 形 等 @9。 雁 片 云 释放 的 组 发 Y 射线 在 大 气 中 输 运 过 程 
与 大 气 中 的 原子 和 分 子 发 生 光 电 效 应 、 康 普 顿 散射 和 电子 对 效应 等 相互 作用 , 这 一 过 程 可 以 
采用 蒙特 卡 罗 方 法 进行 精确 描述 。 

前 人 的 研 完 往往 忽略 碎片 云 的 形状 变化 ,采用 典型 爆 高 典型 当量 对 应 的 半径 和 高 度 经 验 
公 邢 来 描述 组 发 辐射 源 。 对 于 组 发 Y 射 线 在 大 气 中 的 能 量 这 积 ， 则 采用 标准 的 传输 理论 公式 
上 结合 Y 射线 在 大 气 中 的 质量 吸收 系数 来 计算 。 然 而 雁 片 云 半 径 和 高 度 经 验 公 式 不 能 完整 
描述 碎片 云 形状 的 时 空 变化 ， 且 在 不 同 当 量 不 同 爆 局 情况 下 适用 性 较 差 。 男 一 方面 ,标准 传 
输 理 论 仅 考虑 大 气 对 Y 射线 的 袁 减 , 没有 模拟 Y 射线 在 高 空 大 气 中 的 散射 、 吸 收 等 实际 物理 
过 程 ， 急 略 了 7 射线 散射 对 电离 大 气 的 影 啊 。 

为 了 准确 描述 随时 空 演化 的 组 发 Y 射 线 源 在 高 空 大 气 中 的 输 运 过 程 , 本文 建立 了 高 空 核 
爆炸 碎片 运动 演化 的 流体 力学 模型 ,根据 碎片 云 运 动 参数 建 六 了 绥 发 Y 射线 分 层 等 效 体 源 模 
型 ,然后 采用 基于 质量 厚度 抽样 的 蒙特 卡 罗 方 法 模拟 组 及 YY 射线 在 高 空 大 气 中 输 运 并 电离 大 
气 的 过 程 ， 得 到 组 发 y? 射 线 电离 大 气 的 时 空 分布 特 征 。 

1 计算 模型 与 计算 方法 


1.1 碎片 云 流体 力学 模型 
核 爆 炸 碎 片 云 主要 由 核 爆 炸 能 量 驱动 迅速 各 外 膨胀 , 同时 卷 入 高 空 大 气 并 使 之 成 为 碎片 
云 的 一 部 分 。 在 同 外 膨胀 的 同时 ,碎片 云 在 核 爆炸 XX 射线 对 上 下 方 大 气 不 均匀 加 热 产生 的 
压力 梯度 的 作用 下 而 整体 上 升 . 根 据 人 碎片 云 运动 机 理 , 可 以 将 雁 瞩 云 运动 分 解 成 为 膨胀 过 程 、 
上 升 过 程 和 X 射线 加 热 大气 过 程 三 个 方面 。 
雁 片 云 向 外 膨胀 的 过 程 ， 可 以 用 考虑 卷 入 周围 大 气 的 动量 守恒 方程 0 来 描述 : 
dv 


5 全 为 单位 时 间 单位 体积 内 卷 入 周围 大 气 的 质量 源 ，v 为 球 内 某 点 


RR. 
其 中 S=4nR* 守 六 


速度 ，R 为 球体 半径 ， 其 中 M 为 球体 内 的 总 质量 ，A 为 周围 大 气 密 上 度 ，p(7) 为 密度 分 布 函 


数 ，& 为 大 气 分 子 被 卷 入 的 概率 ，p 为 大 气压 强 。 
将 碎片 云 视 为 一 个 整体 , 在 竖 直 方向 上 受到 云 体 上 下 表面 的 压强 差 和 自身 重力 的 综合 作 
用 ， 其 竖 直 方向 的 运动 方程 可 写成 : 
d’Z 
dz 
式 中 ，Z 为 上 升 高 度 ,p 为 大 气压 强 ， 0 为 肆 片 云 面 元 do 的 外 法 线 与 垂直 辣 上 方 辐 的 
夹 朋 : 
碎片 去 上 升 的 主要 驱动 力 来 源 于 X 射线 对 于 高 空 大 气 的 非 均 匀 加 热 , 在 距离 爆 心 距离 R 
处 ，X 射线 经 过 dR 距离 内 空气 吸收 的 能 量 等 于 dR 球 充 内 空气 内 能 的 增加 中 ， 即 : 


4zR2dRep =(e -ee™™)O, (3) 


& 为 气体 内 能 ， 与 大 气温 度 和 p/p 密切 相关 ，0, 为 X 射线 总 能 量 。 联 合式 (3) 和 稀薄 
空气 高 温 状态 方程 03 可 求解 得 到 X 射线 在 大 气 中 沉积 而 引起 的 大 气温 升 。 

根据 式 (1D)~(3) 将 碎片 云 流体 力学 模型 分 为 碎片 云 膨胀 、 碎 片 云 上 升 和 X 射线 能 量 沉积 
三 个 模块 。 首 先 根据 爆 高 、 当 量 、 弹 体 类 型 和 大 气 环境 参数 ， 计 算 X 射线 在 大 气 中 的 能 量 
沉积 及 其 引起 的 大 气温 升 分 布 ， 然 后 根据 初始 的 雁 片 云 参数 ， 数 值 求 解 式 (D) 得 到 碎片 云 下 
一 时 刻 的 径 向 分 布 ， 求 解 式 (2) 得 到 碎片 云 下 一 时 刻 的 上 上升， 循环 迭代 得 到 碎片 云 的 完整 演 
化 过 程 。 

采用 碎片 云 流体 力学 模型 计算 得 到 爆 高 77km 当量 3.8Mt 美国 TEAK 试验 碎片 云 演化 。 
表 1 给 出 了 计算 的 不 同时 刻 的 碎片 云 最 大 水 平 半径 的 中 心 高 度 与 实验 测量 数据 0 的 对 比 。 
计算 结果 与 测量 数据 基本 吻合 , 变化 规律 一 致 。 另 外 计算 得 到 的 碎片 云 由 最 开始 的 球形 逐渐 
演化 为 类 似 于 圆柱 的 倒 梨 形 ， 与 文献 吃 数 据 一致 。 

表 1 TEAK 试验 碎片 云 参数 模拟 与 测量 结果 的 对 比 


Tab.1 The comparison between the measurement and simulation results of debris in TEAK 


= -| pcosGdo — Mg (2) 


时 间 丰 最 大 水 平 半径 /km 中 心 高 度 /km 
观测 数据 “模拟 ”观测 数据 ”模拟 
30 23 27 / 141 
54 50 43 181 193 
96 75 69 245 271 
120 93 84 277 308 
160 99 108 333 357 


187 105 125 360 381 


225 138 149 401 402 
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图 1 模拟 的 不 同时 刻 的 碎片 云 形状 
Fig. 1 The simulated shape of debris at different time 

碎片 云 形状 变化 与 爆 高 当量 密切 相关 。 图 1 给 出 了 模拟 的 爆 高 40km 当量 100kt 和 爆 高 
80km 当量 4Mt 核 爆 炸 不 同时 刻 的 碎片 云 形状 。 当 当量 较 小 爆 高 较 低 时 , 如 图 1 中 爆 高 40km 
当量 100kt， 碎 卢 云 形状 由 球形 演化 为 类 似 于 一 个 上 下 受到 压迫 的 局 椭 球 。 当 当量 较 大 爆 高 
较 高 时 ， 如 图 1 中 爆 高 80km 当量 4Mt， 碎片 云 形状 由 球形 演化 为 典型 的 倒 和 巢 形 。 因 此 在 建 
立 缓 发 ?射线 源 时 必须 考虑 碎片 云 形状 的 影响 。 
1.2 缓 发 Y 射线 源 模型 

核 爆炸 初始 雁 片 云 为 球形 ， 初 始 扩张 为 球形 扩张 。 随 着 雁 方 云 的 扩张 或 上 升 ， 其 形状 随 
之 而 改变 。 由 于 爆 高 与 当量 不 同 ， 售 片 云 或 保持 为 球形 ， 或 演化 为 倒 旬 形 和 椭 球 形 。 为 了 描 
述 人 碎片 云 的 时 空 变化 特征 ,同时 方便 在 蒙特 卡 罗 模 拟 中 抽取 不 同时 刻 的 射线 源 参数 ,本 文 将 
碎片 云 进行 了 分 层 等 效 处 理 〈 如 图 2 所 示 )。 根 据 不 同时 刻 的 雄 片 云 高 度 和 几何 形状 , 将 其 
划分 为 不 同 厚度 近似 圆柱 体 。 将 碎片 云 等 份额 划分 为 n 层 , 每 层 为 在 水 平面 一 定 高 度 上 均匀 
分 布 的 薄 的 圆柱 体 ， 用 中 心 高 度 、 半 径 和 厚度 三 个 参数 撕 述 其 几何 形态 。 在 每 层 进行 均匀 抽 
样 达到 对 不 规则 体 源 抽样 的 目的 .图 2 中 给 出 的 忆 椭 球 性 和 倒 梨 形 等 不 规则 雄 片 云 就 被 等 效 
为 不 同 高 度 的 多 层 缓 及 Y 射线 体 源 。 抽样 过 程 中 首先 根据 各 个 雄 片 云层 所 占 的 份额 , 抽样 出 
Y 射 线 所 在 的 碎片 云层 数 ， 然 后 根据 这 一 层 雁 片 云 的 高 度 、 半 径 和 厚度 ， 然 后 在 这 一 碎片 中 
均匀 抽样 确定 Y 射线 在 该 层 碎 片 云 的 初始 位 置 .采用 分 层 等 效 模型 可 以 很 好 的 描述 碎片 云 的 
形状 变化 特征 ， 分 层 越 多 越 精确 。 


(a) 
扁 椭 球形 碎片 云 


(b) 
倒 梨 形 碎 片 云 


图 2 顶 球 形 与 倒 梨 形 碎片 云 的 缓 发 y 射 线 源 
Fig. 2 The delayed y-ray source from the spheroidicity debris and inverted pear-shaped debris 
裂变 碎片 中 的 放射 性 核 素 平均 每 次 裂变 的 放射 性 随时 间 的 不 断 衰减 ， 其 释放 y 射 线 的 能 
量 释 放 率 可 近似 表示 为 [3]: 


2 (1)=1.407"” (MeV/s: 裂变) (4) 
式 中 +t 的 单位 为 s。 绥 发 y 冉 线 的 能 谱 变 化 较 小 ， 可 以 近似 认为 谱 形 是 一 致 的 [131。 
组 发 y 光 了 于 总 数 根据 式 (4) 与 裂变 当量 来 计算 ， 初 始 方 癌 按照 各 向 同性 抽样 ， 初 始 能 量 按 


照 组 发 Y 射线 能 谱 抽 样 ， 初 始 位置 根 据 Y 射线 源 分 层 等 效 模型 抽样 。 
1. 3 电离 大 气 的 蒙特 卡 罗 模 拟 

组 及 7 射线 在 大 气 中 输 运 并 产生 电离 的 过 程 , 可 采用 蒙特 卡 罗 05 19 方 法 计算 。 由 于 高 空 
大 气 密度 随 高 度 指数 变化 , 采用 传统 的 基于 步 长 的 蒙特 卡 罗 方 法 ,需要 对 高 层 大 气 进行 分 层 
处 理 03， 使 得 模型 较为 复杂 。 本 文采 用 了 基于 质量 厚度 抽样 的 蒙特 卡 罗 方 法 ， 其 核心 思想 
就 是 对 粒子 两 次 碰撞 之 间 的 质量 厚度 进行 抽样 来 蔡 代 传统 的 步 长 抽样 方法 .利用 该 方法 无 需 
对 非 均 匀 连 续 介 质 分 层 处 理 ， 即 可 实现 蒙特 卡 罗 模 拟 。 

假设 第 m+l 次 碰撞 的 位 置 rn+l= rmt+LQm,， 其 中 工 为 两 次 碰撞 之 间 的 输 运 长 度 ， Qnm 为 运 
动 方 辐 。 在 传统 的 蒙特 卡 罗 方 法 中 工 通 过 对 分 布 密度 函数 直接 抽样 得 到 的 。 本 文 则 先 抽样 
得 到 质量 厚度 M〈 即 路 径 上 单位 截面 空气 柱 包 含 空气 的 质量 ) ， 再 由 质量 厚度 M 计算 得 到 
步 长 M。 其 中 质量 厚度 M 的 抽样 公式 51 为: 


Iné 
(CBE) 


(5) 
其 中 ,8 为 (0,1) 区 间 上 的 随机 数 ,5 为 光子 第 m 次 碰撞 后 的 能 量 , 3, (EE,)=2,(E,)/p， 


5%.(E,, ) 为 粒子 与 介质 相互 作用 的 宏观 总 截面 ，p 为 沿路 径 上 的 密度 。 将 质量 厚度 M 转换 


为 步 长 KL， 即 可 完成 非 均 匀 连 续 介 质 中 的 粒子 步 长 抽样 。 

图 3 给 出 了 单 层 雁 片 云 产生 的 组 发 y 射 线 在 大 气 中 输 运 的 几何 模型 。 图 中 ,内 球 为 地 球 ， 
半径 RE 取 6370km。 外 球 是 模拟 区 域 的 上 边界 ， 半 径 Ru 取 7370km。 内 外 球 之 间 是 空气 ， 
空气 简化 为 成 分 固定 的 连续 介质 。 空 气 的 组 成 是 氮 元 素 CN) 占 78%， 氧 元 素 (0) 占 21%。 
空气 的 密度 由 MSISE-00 大 气 模式 0 给 出 。 碎 片 云 半径 Ra， 高 度 hs， 厚 度 Tu， 根 据 分 层 等 
效 模型 取 值 。 采 用 基于 质量 厚度 抽样 的 蒙特 卡 罗 方 法 模拟 跟踪 缓 发 Y 射线 在 高 空 大 气 中 的 输 
运 过程 记 录 源 粒子 在 大 气 中 的 能 量 沉积 分 布 E(r) 。Yy 射 线 在 空气 中 沉积 33eV 的 能 量 便 可 产 
生 一 个 电子 一 离子 对 ， 因 此 由 能 量 沉积 分 布 Er) 即 可 得 到 辐射 在 大 气 中 产生 的 附加 电子 数 


密度 分 布 为 N.(r)= E(r) /33eV 。 


图 3 组 发 y 射 线 在 大 气 中 输 运 的 几何 模型 
Fig. 3 Geometry model of the transport of delayed y-rays in the atmosphere 


对 模拟 中 子 、 光 子 在 大 气 中 输 运 的 蒙特 卡 罗 程 序 MCATNPHS 2 进行 二 次 开发 ， 实 现 了 
基于 雁 片 云 运动 模型 的 组 发 ?射线 电离 大 气 的 数值 计算 功能 。 为 了 验证 程序 ， 采 用 MCNPPR?1 
和 MCATNP 程序 模拟 了 组 发 y 射 线 在 高 空 大 气 中 的 输 运 过 程 。 绥 发 ?射线 均 勺 分布 于 高 度 
100km 半径 10km 厚度 0.1km 的 碎片 云 中 ， 能 谱 取 组 发 y 射 线 能 谱 ， 总 数 取 1, 方 同 为 各 问 同 
性 。 图 4 给 出 了 MCNP 和 MCATNP 程序 计算 的 组 发 y 注 量 和 能 注 量 随 高 度 的 分 布 情 况 。 结 
果 表 明 MCATNP 和 MCNP 程序 的 计算 结果 符合 较 好 。 
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图 4 组 发 Y 射线 注 量 和 能 注 量 随 高 度 的 变化 
Fig. 4 The fluence and energy fluence of delayed y-rays vs. height 


2 计算 结果 与 分 析 

采用 以 上 模型 计算 了 爆 高 40km 当量 4Mt、 爆 高 80km 当量 4Mt、 爆 高 40km 当量 100kt 
和 人 爆 高 80km 当量 100kt 条 件 下 ， 爆 后 不 同时 间 由 碎片 云 产 生 的 组 发 Y 射线 电离 大 气 形成 的 
电子 产生 率 分 布 。 计 算 取 碎片 云 条 件 参 数 与 1.1 节 一 致 ， 取 裂变 一 聚变 比 为 0.5:0.5。 

取 爆 高 80km 当量 4Mt， 取 碎片 云 初始 参数 与 1.1 节 一 致 ， 取 裂变 一 聚变 比 为 0.5:0.5， 
将 碎片 云 按 份额 均匀 分 成 10 层 ， 从 上 到 下 依次 命名 1~10 层 。 计 算得 到 爆 后 Smin 不 同 层 碎 
片 云 释放 的 缓 发 ?射线 电离 大 气 引起 的 电子 产生 率 如 图 5 所 示 。 不 同 层 的 碎片 云 电离 大 气 的 
范围 和 强度 具有 明显 差异 ， 第 1 层 碎片 云 高 度 最 高 ( 约 700km) ， 电 离 大 气 范 围 最 广 ， 而 电 
子 产生 率 峰 值 最 小 ， 第 10 层 碎片 云 高 度 最 低 〈 约 300km) ， 电 离 大 气 范围 最 小 ， 而 电子 产 


生 率 峰 值 最 强 。 总 的 电子 产生 率 等 于 1~10 层 碎片 云 引起 的 电子 产生 率 的 总 和 , 如 图 6 所 示 。 
由 图 可 知 不 同 层 的 碎片 云 电离 大 气 的 强度 与 范围 差异 明显 ， 碎 片 云 形状 对 缓 发 y 射 线 电离 大 
气 范围 和 强度 存在 明显 影响 。 
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5 爆 高 80km 当量 4Mt 爆 后 Smin 不 同 层 碎片 云 引起 的 电子 产生 率 
Fig. $ Electron production rate at 9 minutes caused by the debris in different layer with the explosion height 
of 80km and the equivalent of 4 Mt. 
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6 (b) t=20min 电子 产生 率 10'k 
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6 爆 高 80km 当量 4Mt 爆 后 Smin (a) 和 20min(b) 缓 发 y 射 线 引 起 的 电子 产生 率 
Fig. 6 Electron production rate at $ minutes (a) and 20 minutes (b) caused by delayed y-rays with the 
explosion height of 80km and the equivalent of 4 Mt. 


图 6~ 图 9 分 别 给 出 了 不 同 爆 高 不 同 当 量 情况 下 ， 爆 后 Smin 和 20 分 钟 缓 发 Y 射线 电离 
大 气 形成 的 电子 产生 率 分 布 。 由 图 可 知 ， 电 子 产生 率 左右 对 称 分 布 ， 水 平方 向 上 越 靠 近 爆 心 
电子 产生 率 越 大 ， 垂 直方 向 上 30km 左右 达到 最 大 。 组 发 y 射 线 的 阻止 高 度 约 为 30km， 其 能 
量 主要 沉积 约 在 30~40 公里 的 范围 内 。 这 是 由 于 40km 以 上 大 气 密度 稀薄 ， 其 对 组 发 ?射线 
的 作用 较 弱 。 随 着 高 度 逐 渐 降 低 ， 气 体 随 高 度 降 低 而 指数 增加 ， 对 组 发 射线 的 吸收 急剧 增 
强 ， 从 而 使 得 缓 发 ?射线 主要 沉积 约 在 30~40 公里 的 范围 内 。 


(a) t=5min 电子 产生 率 10"k 


(b) t=20min 电子 产 
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7 爆 高 40km 当量 4Mt 爆 后 Smin (a) 和 20min(b) 缓 发 ?射线 引起 的 电子 产生 率 
Fig. 7 Electron production rate at $ minutes (a) and 20 minutes (b) caused by delayed y-rays with the 
explosion height of 40km and the equivalent of 4 Mt. 
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(b) t=20min 电子 产生 率 10A 


8 爆 高 80km 当量 100kt 爆 后 Smin (a) 和 20min(b) 绥 发 y 射 线 引 起 的 电子 产生 率 
Fig. 8 Electron production rate at 5 minutes (a) and 20 minutes (b) caused by delayed y-rays with the 
explosion height of 80km and the equivalent of 100 kt. 


(a) t=5min 上 


外 子 产生 率 10A 


(b) t=20min 电子 产生 率 10A 


9 爆 高 40km 当量 100kt 爆 后 Smin (a) 和 20min(b) 缓 发 ?射线 引起 的 电子 产生 率 
Fig. 9 Electron production rate at $ minutes (a) and 20 minutes (b) caused by delayed y-rays with the 
explosion height of 40km and the equivalent of 100 kt. 


由 图 可 知 ， 爆 高 40km 和 爆 高 80km 当量 4Mt 核 爆 炸 爆 后 Smin 组 发 y 射线 电离 大 气 的 
电子 产生 率 峰值 分 别 在 107e/s 和 105e/s 以 上 ， 爆 后 20min 电子 产生 率 峰 值 分 别 在 105e/s 和 
105e/s 以 上 。 爆 高 40km 和 爆 高 80km 当量 100kt 爆 后 Smin 组 发 Y 射线 电离 大 气 的 电子 产生 
率 峰 值 分 别 在 105e/s 和 104e/s 以 上 , 爆 后 20min 电子 产生 率 峰 值 分 别 在 104e/s 和 103e/s 以 上 。 
对 比 图 5~ 图 8, 可 以 发 现 : 当量 从 100kt 增加 4Mt 时 , 组 发 Y 射线 电离 大 气 范围 和 强度 增 大 ; 
爆 高 从 40km 增加 到 80km 时 ， 绥 发 Y 射线 电离 大 气 范 围 增 大 ， 但 强度 有 所 降低 。 

图 9 给 出 了 求解 电离 复合 方程 0 得 到 的 爆 高 80km 当量 4Mt 爆 后 Smin 和 20min 的 缓 发 
Y 射线 电离 形成 的 电子 数 密度 分 布 情况 。 电 子 数 密度 峰值 约 在 70km 高 度 附近 ， 这 是 由 于 低 
空 大 气 稠密 ， 电 子 产 生 后 很 快 与 大 气 分 子 发 生 附 着 反应 而 消失 。 根 据 大 气 中 电子 、 原 子 和 分 


子 的 电离 复合 过 程 规律 ， 结 合 电离 复合 方程 求解 的 情况 取 电 子 产 生 率 大 于 10e/s 作为 组 发 
射线 电离 大 气 边界 的 判 据 。 则 有 爆 高 40km 当量 4Mt 爆 后 Smin 组 发 Y 射线 电离 大 气 范 围 接 
近 1200km, 对 应 碎片 云 高 度 约 65km, 半径 约 1Skm, 爆 后 20min 电离 大 气 范围 仍 接近 1000km。 
根据 计算 爆 高 80km 当量 4Mt 焊 后 Smin 组 发 Y 射线 电离 大 气 范 围 可 超过 2000km, 对 应 碎片 
云 高 度 约 400km， 半 径 约 200km， 爆 后 20min 电离 大 气 范 围 仍 超过 2000km。 爆 高 40km 当 
量 100kt 爆 后 Smin 组 发 Y 射线 电离 大 气 范 围 接近 600km， 爆 后 20min 降低 至 500km 以 下 。 
爆 高 80km 当量 100kt 爆 后 Smin 组 发 Y 射线 电离 大 气 范围 超过 1300km,， 爆 后 20min 仍 接近 
1000km。 


(a) t=5min electron density/10" (b) t=20min electron density/10" 
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图 10 爆 高 80km 当量 4Mt 爆 后 Smin (a) 和 20min(b) 的 缓 发 y 射线 电离 形成 的 电子 数 密度 
Fig. 10 Electron density at $ minutes (a) and 20 minutes (b) caused by delayed y-rays with the explosion 
height of 80km and the equivalent of 4 Mt. 


组 发 y 射 线 电离 大 气 强 度 与 范围 随 爆 高 、 当 量 和 时 间 变 化 复杂 。 组 发 y 射 线 源 与 碎片 云 参 
数 密切 相关 。 一 方面 碎片 云 受到 X 射线 对 爆 心 下 方 低层 大 气 加 热 后 的 浮力 推动 上 升 ， 使 得 
绥 发 y 射 线 源 不 断 上 升 ， 进 而 使 得 电离 大 气 范围 增加 ， 电 离 大 气 强 度 降 低 。 男 一 方面 绥 发 
射线 源 强 度 随 着 时 间 指 数 衰 减 ， 但 碎片 半径 随时 间 增 加 。 

当 爆 高 相同 当量 不 同时 ， 核 爆炸 缓 发 y 射 线 源 强度 与 当量 呈 线 性 关系 。 随 着 当量 增加 ， 
射线 总 数 增加 ， 碎 片 云 高 度 和 水 平 半径 增 大 ， 组 发 y 射 线 电 离 大 气 范 围 增加 ， 强 度 峰 值 一 般 
也 有 所 增加 。 对 比 图 5~ 图 8， 当 量 由 100kt 增加 至 4Mt 时 ， 组 发 y 射 线 电离 大 气 范 围 由 几 百 
公里 增加 至 千 公 里 以 上 ， 强 度 峰 值 增加 1~2 个 数量 级 。 当 当量 相同 爆 高 增加 时 ， 碎 片 云 高 度 
和 水 平 半 径 增 大 ， 组 发 ?射线 电离 大 气 范 围 增 加 ， 但 射线 总 数 保持 不 变 ， 所 以 幅 值 会 有 所 降 
低 。 对 比 图 5~ 图 8， 爆 高 由 40km 增加 至 80km 时 ， 缓 发 y 射 线 电 离 大 气 范 围 增 加 约 一 倍 ， 强 
度 峰 值 降低 1~2 个 数量 级 。 

3 总 结 

为 了 准确 描述 随时 空 演化 的 组 发 Y 射线 源 电离 大 气 过 程 , 本 文 建立 了 高 空 核 爆炸 碎片 运 
动 演化 的 流体 力学 模型 ， 然 后 根据 碎片 云 运动 参数 建立 组 发 Y 射线 分 层 等 效 体 源 模型 , 通过 
分 层 抽 样 得 到 组 发 Y 射线 源 的 初始 位 置 、 能 量 , 采用 基于 质量 厚度 抽样 的 蒙特 卡 罗 模 拟 方法 
计算 得 到 组 发 Y 射线 电离 大 气 的 电子 产生 率 。 计 算得 到 不 同 爆 高 不 同 当量 情况 下 爆 后 不 同时 
刻 组 发 yY 射线 电离 大 气 形成 的 电子 产生 率 分 布 情况 。 碎 片 云 形 状 对 组 发 Y 射线 电离 大 气 范围 
和 强度 存在 明显 影响 。 研 究 得 到 了 组 发 Y 射线 电离 大 气 强度 与 范围 随 爆 高 和 当量 的 变化 特征 
规律 。 
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